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O '■ Abstract 

r~| ■ Inflationary correlation functions are potentially affected by infrared divergences. For example, 

Qh. the two-point correlator of curvature perturbation at momentum k receives corrections ~ ln(kL), 

where L is the size of the region in which the measurement is performed. We define infrared-safe 
correlation functions which have no sensitivity to the size L of the box used for the observa- 
tion. The conventional correlators with their familiar log-enhanced corrections (both from scalar 
and tensor long-wavelength modes) are easily recovered from our IR-safe correlation functions, 
■^j- ' Among other examples, we illustrate this by calculating the corrections to the non-Gaussianity 

parameter /tvl coming from long-wavelength tensor modes. In our approach, the IR corrections 
automatically emerge in a resummed, all-orders form. For the scalar corrections, the resulting 
' all-orders expression can be evaluated explicitly. 

p : 

^ ; 1 Introduction 

,— 1 \ Infrared divergences associated with the inflationary power spectrum are a long-standing issue [IHS] which 
' has more recently received a lot of attention following [3 IB]- Our focus is on divergences which directly 
affect the power spectrum. According to [9] such divergences are cut off by the size L of the observed 
5_i , patch of the late universe. We are going to develop and generalize the analyses of [lQl[TTJ (see also [12]). 
This approach emphasizes the way in which long-wavelength perturbations do (or do not) influence locally 
measured inflationary spectra. It may be related to earlier proposals of |13]I14] and is clearly in line with at 
least part of the subsequent discussion in [15] . 

Infrared divergences explicitly arise in loop corrections to inflationary observables (using e.g. the in-in 
formalism) or through the nonlinear dependence of curvature perturbation on fluctuations of an underlying 
scalar field (e.g. in the 5N approach). The vast amount of literature on the subject (see e.g. |16H4U| ) has 
recently been reviewed in [41] . More or less by definition, IR divergences are due to modes which have 
a much longer wavelength than the characteristic scale of the problem under consideration. Focussing on 
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correlation functions, it is then clear that such IR modes left the horizon earlier than the modes which are 
directly accessed via the correlation function at a given scale. Hence a generic effect of these IR modes is a 
modification of the background in which other modes propagate and are eventually observed. 

Much debate has been raised on how to interpret such IR divergences. At least in the case of single-field 
slow-roll inflation, it is now widely accepted that log-divergent integrals over soft modes have to be cut off 
at a k m i n ~ 1/L, where L is the typical size of the 'box' in which the observer measures some correlator 
at momentum k 3> 1/L [9] (see also [42H45j ) . It is thus more appropriate to talk about a log-enhancement 
rather than a log-divergence. The suggestion that long-wavelength modes can be absorbed in the background 
and hence do not affect the correlator at momentum k has been put forward long ago |46H51j . In this sense, 
the absence of IR divergences in situations where L is not much larger than k may have been apparent to 
many authors even before Lyth's paper [U] of 2007. 

We note that backreaction of long-wavelength modes in quasi-de Sitter space-times has been considered 
also in other contexts |48H60| . most notably in attempts to compensate the cosmological constant or to 
explain the current accelerated expansion of the Universe. We have nothing to say concerning these topics 
and refer the reader to [H] for a more extensive compilation of the relevant literature. 

In this paper, we analyse IR effects associated with the backreaction of long-wavelength scalar and tensor 
modes in inflationary backgrounds, in the spirit of [1Q|. I1 Ij . This approach is related to the consistency 
relations [6TH63] . Developing and generalizing a suggestion made in [11] (see also [E]), we propose an 
IR-safe definition of correlation functions involving curvature fluctuations. In doing so, we remove any 
sensitivity to modes that have a much longer wavelength than the scale at which the correlator is probed. 
The essential idea is to make use of the proper invariant distance on the reheating surface, where the 
curvature perturbation is evaluated. Considering two points on this surface, the dependence of the (physical) 
invariant distance on the coordinate vector, corresponding to the separation of the two points, is affected by 
long-wavelength contributions from geometrical quantities, namely the curvature and tensor perturbations. 
The misidentification of the distance due to long-wavelength modes is precisely the origin of IR effects. 
Consequently, by using the proper invariant distance, it is possible to construct n-point functions for the 
curvature perturbation that are free from the effect of long-wavelength modes and, hence, free from IR 
divergences associated with these contributions. 

We show how to relate n-point functions, calculated in terms of the invariant distance, to the conven- 
tionally defined n-point functions. This allows us to provide closed expressions for the latter that manifestly 
exhibit the dependence on long-wavelength modes. As a consequence, in our approach the IR corrections to 
n-point functions automatically emerge in a resummed, all-orders form. When expanded at leading order 
in terms of long-wavelength modes, we recover the familiar log-enhanced, IR sensitive contributions. We 
apply our approach to the analysis of the two- and three-point functions for the curvature perturbation in 
slow-roll, single field inflation. The leading IR corrections to the power spectrum appear as log-enhanced 
contributions, multiplied by the power spectrum and second order slow-roll parameters. Furthermore, our 
resummed, all-orders expression allows us to evaluate IR corrections in a non-perturbative way by using 
statistical properties of the integrated long-wavelength fluctuation. We apply this framework to specific in- 
flationary set-ups, obtaining a complete expression that includes all contributions of scalar long- wavelength 
modes to the power spectrum. Regarding the bispectrum, we derive the complete expression for long- 
wavelength scalar and tensor contributions to /nl- We then expand the result at leading order in slow-roll, 
showing that tensor modes dominate the slow-roll expansion and provide the leading log-enhanced contri- 
butions to non-Gaussianity. Contrary to the power spectrum, we find that the leading order correction to 
/nl is suppressed only by first order slow-roll parameters. 

We also show that, in all cases where the <5iV-formalism is applicable, our results can be equivalently 
obtained in terms of a suitable generalization of the <5iV-formalism, extending the discussion of [11]. In the 
present work, we include the effects of graviton long- wavelength modes, and we explain how to calculate IR 
contributions to arbitrary n-point functions involving curvature perturbations. 
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Log-enhanced contributions to inflationary observables, both in the in-in formalism and in the 6N- 
formalism, have received much attention over the past few years. In the case of the <5iV-formalism, they 
have been associated with infrared divergences of the so called C-loops, and have been calculated mostly 
in terms of a diagrammatic expansion [43 ^641455] . Although our approach is related, it is conceptually and 
technically different. We derive IR contributions directly from geometrical quantities. These corrections 
appear automatically in a resummed, all-orders form and do not need any diagrammatic expansion. Fol- 
lowing arguments given in |10| . the presented derivation of IR effects from the geometry of the reheating 
surface matches IR contributions due to quantum loop effects of long-wavelength modes calculated a la 
Weinberg [7JIE], although a direct comparison is beyond the scope of the present paper. 

The paper is organized as follows: In sec. [5J we present the definition of IR-safe correlation functions and 
we show the emergence of IR corrections in a resummed, all-orders form by relating conventional correlation 
functions to their IR-safe equivalents. Furthermore, we show for scalars how this expression can be evaluated 
explicitly. In sec. [3j we give an alternative approach in single field, slow-roll inflation in the language of the 
5 N- formalism. In sec. 0]and sec. [5J we apply our results to the power spectrum and to the non-Gaussianity 
parameter /nl, respectively. We draw our conclusions in sec. El 



2 Geometry of the reheating surface 

In this section we provide a physical interpretation for the appearance of log-enhanced correction to infla- 
tionary correlation functions, developing and generalizing [10[lll|. We start from the assumption that the 
reheating surface (or any other surface of constant energy density after the end of inflation), viewed as a 
metric manifold, represents in principle a physical observable. Neglecting vector modes, the metric of this 
surface can be written as 

dsl = e 2 ^ s) (e^) dx i dx j . (1) 

We choose a gauge where the symmetric matrix 7 is traceless and dijij = 0. Of course, £ is accessible only 
indirectly, e.g. via 5T/T of the CMB radiation, but for the present paper we simply assume that this does 
not limit its observability. Considering, for example, single field slow-roll inflation, the expressions for £ and 
jij (in momentum space) read [67] 

N^(q)H(q) ^ H(q) s s 

V^Q s=+,x vQ 

These quantities are conserved on superhorizon scales |6 H I68 |, [69] . In the equations above, and are 
normalized Gaussian random variables and s is the helicity index for gravitational waves. The polarization 
tensors ef • are chosen to satisfy the transversality and tracelessness conditions, as well as an orthogonality 
relational Furthermore, N^q) = V/(dV/d<j>) and H(q) = y/V((f>)/3 with both quantities evaluated at the 
time of horizon exit of the mode q. 

Although, for definiteness, we focus on slow-roll inflation, the particular expressions for £ and 7^ given 
above are not essential for the formalism presented in this section. Consequently, our arguments are largely 
independent of the specific inflationary set-up under consideration. Important consequences for the n-point 
functions can be derived just from the geometry of the reheating surface specified above. Focussing on 
corrections to the curvature perturbation £, we start by discussing the power spectrum and then generalize 
to spectra of n-point functions, for arbitrary n. 



4 We use conventions such that the Fourier transform reads £(x) = J e lkx C,{k) . The Gaussian random variables and 

b~ have zero mean and variance (ag ap) = (27r) 3 &^{k + p) and (6| bp ) = {2tt) 3 6^ 3 '(k+p) S ss . The polarization tensor for 
gravitational waves satisfies efj(fc) = = kiefj(k) and the orthogonality relation £V efj(k) efj(-k) = 28 ss i . 
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2.1 The power spectrum 



Using its definition, (£j Cp) = (2ir) 3 5^{k + p) 2vr 2 T(;(k)/k 3 , the power spectrum can be written as the 
Fourier transform of the correlation function in real space: 




Since we want to interpret this formula as a practical prescription for the measurement of the power spectrum, 
we do not view (• • • ) as an abstract ensemble average. Instead, the averaging is over pairs of points separated 
by a coordinate- vector y within a certain part of the reheating surface. In other words, we are averaging over 
the location x of such pairs. This prescription clearly relies on a certain parameterization of the reheating 
surface and is hence gauge dependent. Nevertheless, given the gauge choice made earlier, the resulting V^(k) 
is a well-defined observable. 

An observer is not able to probe the whole inflationary region. We assume that the observable patch 
is a box of volume L 3 to which the x-averaging is restricted. While this may not be immediately apparent 
from ([3D, the power spectrum measured by a given observer depends on the box-size L. Qualitatively, this 
can be seen as follows: 

Focus on a certain momentum k. Due to the Fourier transform, the power spectrum at this k is 
determined by the behavior of the ( ((x) C(x+y)} as a function of y in the region y ~ 1/k (here y = \J 5ij y % yi 
is the length of y). However, at different x the same value of y may correspond to different physical (invariant) 
distances between points x and x + y at which £(x) and £(x + y) are evaluated. The reason for this is the 
long-wavelength background 



which varies significantly as x varies over a box of (sufficiently large) size L. Indeed, the physical distance 
between the points x and x + y appearing in the average is given by z 2 = e 2 ^ (e^- y l yi ' . Moreover, this 
mismatch between y and the true distance z grows with L. This effect is at least one of the origins of the 
familiar IR-problems of inflationary correlation functions. At leading order, IR-problems originate precisely 
from this effect. 

To be more precise, we somewhat jump ahead and note that (£ 2 ) ~ (N^H) hi(kL), with a similar 
formula holding for 7. In other words, the expectation value of Q 2 grows logarithmically with L because 
of the summation over modes between 1/L and k involved in its definition. Thus, the effect of these 
backgrounds on ^-correlators at the scale k can become large if the logarithm overcomes the suppression 
by the tree-level power spectrum (N^H) 2 . Such a potentially large effect can come only from the factors 
e 2 ^ and relating the coordinate distance y and the invariant distance z, as explained above. If we are 
able to remove this effect from the definition of the power spectrum, then we have removed all IR effects at 
the leading-logarithmic order. By this we mean all corrections involving as many powers of \n(kL) as of the 
suppression factor (N^H) 2 , at leading order in slow-roll. 

To avoid this (leading-logarithmic) L-dependence (or IR-sensitivity) , we propose to use the invariant 
distance z for the definition of the curvature correlator [1 1|. I15j . The background contains, by its very 
definition, only modes much longer than the relevant scales y ~ 1/k. Hence the background is smooth at 
the scale y. Its presence corresponds to a (constant) coordinate transformation y — >■ z: 






L- 1 <q<s:k 



L~ 1 <q<^k 




J 



(5) 



The invariant distance z = yj 5ij z l zi represents the physical separation of the points x and x + y , indepen- 
dently of the location x and the background in its surroundings. 
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Thus, the correlator (((x) £(x + e~^> e -7 ^/ 2 z) ) involves an average over pairs of points that are 
separated by a certain invariant distance z. The z-dependence of this correlator is then a background- 
independent object. To make this even more apparent, we spell out the exact prescription for obtaining this 
correlator: The basic step consists in picking a pair of points from the reheating surface which are separated 
by an invariant distance z and multiplying the corresponding values of £. This in itself is not background 
independent since the background can shift ( by a constant. However, once we restrict our interest to 
the z-dependence of this product of C-values, any such constant drops out. Hence, the z-dependence of 
( C(x + e - ^) e - ^)/ 2 z) ) is indeed an IR-safe quantity: While the average is in practice over a certain 
region of size L, the expectation value is independent of where we are in this region. It can therefore not 
depend on the size L of the underlying region. To say it yet in another way: Single-field inflation ends in 
the same way in every part of the universe and hence a correlator, defined in a purely local manner, can not 
depend on the size of the region from which the sample of pairs of points is chosen. 

Consequently, we can define an IR-safe power spectrum, that we denote , by 

3 



Vf\k) = ^ J d 3 z ( ({x) ((x + e-tW e" 7 ^)/ 2 z)^ . (6) 

This Fourier transform at scale k is only sensitive to the z-dependence of the correlator in the region z ~ 1 jk. 
It is hence IR-safe by the arguments given above. 

The expression for the original IR-sensitive power spectrum given in ^ follows from comparing 
eq. ([3]) and eq. ([6]) . Starting from eq. ([3]) , we express the vector y in terms of the vector z. Notice that this 
also affects the argument of the exponential. Then, we perform a coordinate transformation d 3 y — > d 3 z in 
order to bring the integral in a form similar to eq. ([6]). As a final step, we can express the result in terms of 
the IR-safe power spectrum evaluated at e^"e~^™ 2 i . A detailed calculation can be found in Appendix 
[Al The result reads 

P c (k) = 



1 -3/2 \ 

(e-^)).>% 7f> ( e -«*> e -tf*>/*£) \ . (7) 



The vector k is a unit-vector in fc-direction and the average is performed over the background quantities £(x) 
and jij(x). Neglecting tensor fluctuations in the equation above, we recover our result [11] for corrections to 

the power spectrum due to scalar fluctuations, i.e. V((k) = (V^ (ke~^) ). Let us point out the presence of a 
prefactor, containing only tensor fluctuations, in eq. 0. It is originating from the coordinate transformation 
d 3 y — > d 3 z in the comparison of the two spectra. While scalar fluctuations receive a contribution from this 
transformation, tensor fluctuations do not, due to the fact that det e 7 = 1. Expanding to leading non-trivial 
order in the background yields 

nW=(l-^(trf)^ + i(?)^)<W , (8) 

in agreement with [10] (see also Section H]). Here, we used the zero mean condition (((x)) = = (^(x)), 
which can always be realized by a rescaling of coordinates. In principle, one may choose coordinates where 
this is not the case. But it is rather natural to assume that an observer would specify coordinates in such 
a way that his observable patch is not affected by a constant background shift. In the particular case of 
slow-roll inflation, both corrections in eq. ([8|) are of the same order. While, according to the scalar-to-tensor 
ratio, tr7 2 is more slow-roll suppressed than Q 2 , it appears with only one derivative in In A;. Hence, tensor 
corrections are as important as scalar corrections in slow-roll inflation. 

The remaining task is to average the background quantities, given in eq. In principle, we have 

to average (,(x) and j(x) over the large observed region of box-size L. However, this is equivalent to an 
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ensemble average of £(0) and 7(0) with IR cut-off L. Thus, in single-field, slow-roll inflation, we are dealing 
with sums of Gaussian random variables a^, respectively bt, 

d 3 q ^ f d 3 q N^q) H(q) 



{^r J ( 27r ) v2o 3 

1/L 1/L 



k k 



* - / ■»<*> -!§?^ H H ^ * • (10) 

1/L 1/L S_+ ' X 

While their averages are vanishing, (£) = = (jij), one finds a scale-dependent result for the two-point 
functions. For instance under the assumption of a scale-invariant power spectrum, they obey a logarithmic 
scale-dependence 

(C 2 ) = (^pfHkL) (tr7 a > = (7«7«>=8^) ln(fcL) . (11) 

Neglecting tensor fluctuations for the moment, the background Q is a sum of Gaussian random variables a^. 
Therefore, Q itself is a Gaussian random variable, with distribution^ 

F[C]dC= — L= exp f-iL) dC (12) 



where the width is 



^ V 2CJ C, 



d 3 g N 2 {q)H 2 {q) 



l/L 

Note that we do not assume a scale-invariant behavior of the power spectrum in this expression. Typically, 
the n-point functions we are interested in can be expressed as ( /( £(x) ) ), for some function f. As usual for 
Gaussian variables, this may be expressed in terms of an integral over a Gaussian probability distribution 

(/(C)) = -iL= I d( exp f -il ] /(C) • (14) 

As an example, the power spectrum is given by 

V ( (k) = -jL= [ d( exp vf(ke^) . (15) 

^a 2 J V 2a Cj 

Consequently, the question of convergence of fluctuations due to long-wavelength modes reduces to con- 
vergence properties of this single integral. The usual series expansion can be recovered by expanding the 
function in the logarithm of the scale k. This yields 

~ ,72n\ d 2n Vf\k) 



n=0 

( ( 2n > = ~jL= I d( C 2n exp ( ) = (2n - 1)!! (a 2 ) n , (17) 



= / dC( 2n exp(-^ 



3 This is related to the stochastic approach [TB] of Starobinsky. 
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where n!! denotes the double factorial. This is in agreement with [10]. We emphasize, however, that a 
breakdown of convergence of the series does not necessarily mean a breakdown of convergence of the integral 
in eq. (|15p . We return to this point in sect. [H Notice also that in eqs. ([9]) and ()10p we have neglected the 
intrinsic non-Gaussianity of curvature and tensor perturbations. Such intrinsic non-Gaussianity is present 
at sub-leading order in slow-roll. However, at every log-order, there is a term consisting solely of Gaussian 
contributions. Relative to this term, contributions with intrinsic non-Gaussian parts are suppressed by 
slow-roll parameters and the Hubble scale with no additional log-enhancement. Therefore, neglecting the 
non-Gaussian contribution is justified in our leading-log analysis of IR-corrections. 

Attention has to be paid to the fact that inflation has ended at some point. Hence, there exists a value 
kmax corresponding to modes that have never left the horizon. The observer measuring V^{k) for some fixed 
k will have to exclude regions where ke~^ > k ma _ x from his averaging procedure. Technically, this implies a 
lower bound for the (- integral, given by ( mm = — hi(/c max /fc). 




Vdk) = ~^= 2 J dC exp ( ] 7>w(fce-') (18) 

The Gaussian function in eq. (|18p gives a non-negligible contribution only in a limited range around zero. 
This range is of the order of o~^. For large L (implying large a^) and for k sufficiently close to fc max , the lower 
bound Cmin enters this range. Hence, in such cases, the lower bound implies the subtraction of a significant 
contribution from the integral. We finally note that the existence of /c max and Cmin are related to potential 
convergence problems of the series expansion in eq. (|16p . This is apparent since the slow-roll conditions, 
which are responsible for the smallness of derivatives of V^(k), break down near fc max . 

Including tensor modes is in principle straightforward, but complicated by the matrix structure of 7 and 
the different independent polarizations involved. In order not to overburden formulae, we set the scalar 
background ( to zero in what follows. The complete power spectrum can then be expressed as 



-3/2 



V C {k) = ( (e-%fc^ "~vf ] ( (e-%fci% *k) ) , (19) 



where k is the unit vector parallel to k. We also introduce the notation 



n 



) lm kkm ■ (20) 



Note that each entry of the matrix 7^, being a sum of Gaussian random variables, is a Gaussian random 
variable. However, the various entries in the matrix are not statistically independent: this implies that it is 
not obvious how to calculate the statistical distribution of the entries of the exponential of (—7), that enters 
in the definition of n. Having this distribution, that we denote with P[n], it is straightforward to provide 
an integral representation for the power spectrum subject to tensor background modes: 



P c (fe) = J dnF[n] n - i(jfe)7^ 0) (n? k) . (21) 
It is clear that, at least numerically, P [n] can be determined and the integral can be calculated. 
2.2 Higher correlation functions 

To discuss n-point functions, we could try to generalize the 'almost scale-invariant' spectrum of eq. ([3]) by 
writing 

r n (k u ...X) = ( (£2) J <i 3 yi • • • d 3 |/» e-«(^' + '- +fe; ^) ((f) ((£ + fi) . . . ((f + i) |) . (22) 
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However, it is not clear which particular combination of k\ . . . k n one should use to define k in the prefactor 
k 3n . This is not irrelevant since factors e 1 ! 2 will get tangled up in this prefactor. Hence, we choose to write 
the general formula for the higher-order analogue of the conventional spectrum PAk) = 2n 2 V^(k)/k s . In 
doing so, prefactors will arise from the scaling of the d?y a . Since the determinant of the tensor contribution 
is one, this scaling consist exclusively of which only depend on the overall scale. Given these preliminaries, 
the generalization of our formalism is completely straightforward and the IR-safe spectrum is defined as 

pW(k u . . . , k n ) = I J d 3 z, . . . d 3 z n e -Kk^i+-+knz n ) c(f) . . . ((x + y n ) \ , (23) 

where 

Va = y a (z a , C, 7) = e" C V^ 2 i* a . (24) 

This means that we measure the correlation function in terms of n invariant distances, characterized by a set 
of vectors z a , a G [1, . . . , n]. Hence, the z a -dependence of the corresponding n-point function is independent 
of background quantities and, therefore, IR-safe. Consequently, its Fourier transform, i.e. the spectrum P^\ 
is the desired IR-safe spectrum. A straightforward generalization of the previous calculation for the power 
spectrum provides the following result 

P n (h,...,k n ) = ( e - 3n ^(°He-^ 2 fc ll ... ) e-^ 2 4)) ■ (25) 

As already stressed above, the prefactor e~ 3n ^ originates from the naive scaling P^ ~ /c _3n . 

The log-enhancement-effects of higher correlation functions specified by eq. (|25p can be directly applied 
to observables measuring non-Gaussianity, like /nl> as we are going to discuss in section [5j 



3 An alternative approach within slow-roll inflation 

In the previous section we discussed a systematic way to define IR-safe n-point functions. We have explained 
how to straightforwardly obtain, from these IR-safe quantities, the corresponding IR-sensitive objects. In 
this section, we present an alternative point of view: working only in momentum space, we will directly 
calculate the all-orders IR-enhancement of the conventional power spectrum. To be more specific, we will 
compute the curvature perturbation £, by implementing a suitable extension of the JA^-formalism, in such a 
way as to include the effects of long-wavelength modes. The results coincide with what we obtained in the 
previous section, in all cases in which ^A^-formalism is applicable. So, in these cases, the two methods are 
equivalent. 

We focus on a single, slowly rolling scalar field 4> (the extension to multiple fields is outlined in Appendix 
IB"]) . We assume the underlying metric to be of the form 

ds 2 = -dt 2 + a 2 (t) g ij dx i dx j . (26) 

Throughout this section, we are interested in quantities evaluated at some wave vector k. To analyze 
contributions from the background to these quantities, we find it technically convenient to separate the 
fluctuations into modes characterized by momenta larger, and smaller, than k. For modes q with q = \q\ <^ k, 
we will work in a gauge with 5^<p = 0. While for q around k and larger, we adopt a gauge with vanishing 
scalar metric fluctuations. The advantage of this splitting, and of these different gauge choices, is that the 
contribution from long-wavelength modes is contained in geometrical quantities and, therefore, contained in 
the 3-metric 

= ^ (e% ■ (27) 
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Here, the scalar and tensor background, £ and 7™, are denned as before. It would be interesting to understand 
whether the above construction can be done in a gauge invariant manner. On the other hand, let us stress 
that we proceed in this way only for technical convenience. One could also work with a gauge characterized 
by vanishing scalar metric fluctuations for all q. With this choice, however, the scalar background from 
long- wavelength 5^<j) would affect the scalar field value at the time of horizon exit (see [11] for a treatment of 
background modes of the scalar field <f> in this latter gauge choice). In contrast, the tensor background would 
still enter via the 3-metric. Therefore, the inclusion of tensor background modes within the <5iV-formalism 
requires a treatment as outlined in this section, contrary to the scalar background which might be calculated 
by different techniques. 

The appearance of background contributions in the 3-metric eq. (|27p has important consequences for 
the physical length scale associated with the wave vector k, i.e. on the physical wavelength. Due to the 
deviation of g~ij from flatness, this scale is not the inverse of k = ^/kikjdij, but is instead given by 1/k' with 

k' 2 = e- 2f (e-% hkj . (28) 

Hence, the physical scale depends on the original vector k and on the background quantities ( and jij. This 
dependence on background quantities leads to a shift in the time of horizon exit for a given scalar mode of 
momentum k, from t k to ty. That is, since the time of horizon exit is defined by the relation k = a(tk)H{tk) 
(a being the scale factor), at first order in slow-roll we have the relation 

dt k = — d\nk . (29) 
H 

For small time variations, and at leading order in slow-roll, we can integrate the previous equation and find 

k! 

H (t k > - t k ) = hi — = -C - A with (30) 
k 

1 1 \ t t fl 



A = l^-Jij - jliHijj hkj + I hkj) + 0{f) . (31) 

Here, k represents a unit vector in /c-direction. The quantity A collects the leading order contributions from 
the long- wavelength tensor modes, obtained from expanding the exponential in eq. ()28|) . 

The form of the background metric affects the dynamics of first order, massless scalar fluctuations. In 
momentum space, the equation of motion for the scalar perturbations reads 

k' 2 

(^)" + 3ff(^)- + ^^ = , (32) 

where dots denote derivatives with respect to time. Note that the effect of background quantities enters via 
the Laplacian which leads to the k' 2 instead of k 2 in the third term on the left-hand side. The solution for 
the fluctuation S^cp results in 

/ H(k') 
6#f> = S^(k', g^) = j— a % (33) 

(det* fy) (2fc'3)2 

in a superhorizon regime. In our notation H{k') indicates that this quantity is evaluated at time of horizon 
exit of the scale k', instead of k, in order to take into account the shift due to long-wavelength contributions. 
The normalization of S^cj), det 1//4 gij in the denominator, is obtained when imposing the usual commutation 
relations between the quantized scalar fluctuation and its momentum conjugate (see, for example, |70j ) . 
Another way to understand it is the following: the normalization of 5^(j) is set by requiring that in the limit 
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of short distances y, we recover the singularity of the scalar field in Minkowski space for ( 5(f>(x) 5(f>(x + y))- 
On these distances, the background quantities Q and 7jj are constant and, hence, need to be absorbed in a 
redefinition of space variables in order to bring the metric in Minkowski-form. This redefinition is responsible 
for the factor det 1//4 gjj, appearing in the normalization of 5j:4>. 

Since jij is traceless, det g~ij = exp(6£) and eq. ([33|) can be rewritten as 



H(k') 



V2 

in! (k) 



hkj (e-^ 
H(k') 



{2k 3 ) 



I H 



(34) 



where we define the function m(k) that depends on a unit vector k along the direction of k: 



m(k) 



(35) 



In equation (j34[) . the dependence on background quantities is limited to the overall function m(k) (that 
depends only on the tensor background, see eq. (|35|) ) and to the 'time' argument k' of the Hubble parameter. 



Starting from scalar fluctuations and by using <5./V-formalism [42 | f51J |l71rf73j . we can express the curvature 
fluctuation £ at superhorizon scales on a constant energy density slice, that we take to be the reheating 
surface, in terms of 5(j). The curvature perturbation £ ? is related to the time integral of the local expansion 
parameter, providing the number of e-foldings, from an initial hypersurface (that we take at time of horizon 
exit for the mode k) to the final hypersurface of constant energy density. In single field inflation, we have 



C = JV [<f> + 84>) - (N) 



(36) 



where (N) is the spatial average of the first term on the right-hand side. The quantity 4> + 5(J) corresponds to 
the homogeneous value for the scalar field plus its perturbation built, as above, on a space-time geometry that 
includes the contributions of long- wavelength modes. The previous schematic expression can be expanded 
in the scalar fluctuations, and gives in momentum space H 



= N (j> (k')5 j; (f>(k , ,g ij ) + 



(37) 



Notice that functions on the right-hand-side are evaluated at time of horizon exit of the mode k, which is 
sensitive to the change in the background geometry due to long-wavelength modes. That is, their argument 
is k' instead of k. As in section [21 the function Nj, = dN/d(j) is given by = V/(dV/d(j)). The remaining 
terms in the 5N expansion, understood in the dots of eq. (|37|) . are slow-roll suppressed with respect to 
the first one. Using the results obtained earlier, we get for ^ an expression in terms of Gaussian random 
variables as follows 



(2k 3 ) 



(38) 



The dependence on long-wavelength background quantities is contained in the overall factor between squared 
parenthesis. Eq. (|38p . possibly including higher order terms in the <5iV-expansion, is all what we need to 
straightforwardly compute inflationary observables, associated to n-point functions of curvature perturba- 
tions, including the effects of long-wavelength modes. Eq. (|38|) can be regarded as an extension of 5N- 
formalism. It includes the contributions of long-wavelength scalar and tensor modes in the expression for 
the curvature perturbation £. 



For the purposes of this work, we can truncate the 5N expansion to the first, leading order term in slow-roll. Including 
higher order terms is straightforward, as we discuss in Appendix iBl 
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As an application of eq. (|38p . we rederive the expression for log-enhanced contributions to the power 
spectrum. We start with the two-point function of the curvature perturbation 

(CkCp) = — (m l 2(k)m l 2(p)N4k')H(k')N4p')H(p')a l :a f ) 

= ^»M (ra(M(W> , (39) 

where for passing from first to second line, we used Wick's theorem and contracted the Gaussian variables 
and ap. Indeed, is only allowed to contract with ap, since any other quantity depends on modes with 
momenta much smaller than k. One obtains the following expression for the power spectrum: 

n« = ^(m(fc)JVj(fc')i? 2 (fc')> • (40) 

The argument of the average on the right-hand side depends on long-wavelength scalar and tensor contri- 
butions, which, as shown in eq. have non- vanishing two-point functions. In absence of contributions of 
long-wavelength modes, eq. ([30]) provides the following tree-level result 

Vf\k) = ^Nl(k)H 2 (k) , (41) 

which also coincides with the definition of the IR-safe power spectrum provided in section [2j Notice that 
the dependence on the scale k in vjP occurs only through the dependence on the time of horizon exit of 
the right-hand side. Using this fact, eq. (|40|) can be rewritten as 

V C (k) = (m(k)vf(k')) . (42) 

Recall that k represents the unit vector along the direction of k, while k' in the previous expression is 
associated to k via eq. (|28p. Using these formulae, eq. (142 p can be rewritten as 

V ( (k) = / [(e-% kik 3 ] ~ 3/2 (e-^l 2 k) \ . (43) 

Not surprisingly, this corresponds exactly to equation (J7J), obtained with the method of section [2J 



4 Two-point function and the power spectrum 



In this section, we analyze the log-enhanced corrections due to scalar and tensor long-wavelength modes to 
the power spectrum of curvature perturbation. Using the results from sec. [2] or sec. El the power spectrum 
is given by the formula 



V C (k) = {m{k)vf ] {k')) 



-3/2 



,(0) 

C 



<e~^ 2 k 



(44) 
(45) 



As explained in section [2j this implies that we can deal with scalar perturbations to all orders, resumming 
the complete series, in case an exact expression for the tree level power spectrum is known. We will return 
to this important topic at the end of this section; for the moment we focus on calculating, in full generality, 
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the leading log-enhanced contributions to the power spectrum. In order to do so, it is sufficient to expand 
eq. (|45|) in £ and 7. The following equations are useful for this purpose 



m(k) 
Ink' 



In k 
1 



c 



k^ kj 
A 



-3/2 



1 + kikj 



3 r r 15 



-it 1 - - 11 1 
Jijkikj ~^)u1lj kikj + — 



(46) 
(47) 
(48) 



Here, k! denotes the Euclidean length of the vector e 7 / 2 A?. We will also make use of the identity (see 
also GDI): 



( lij llr, 



1 

30 



( tr 7 2 ) [3 (SuSjm + 8i m 6ji) - 2 5ijS lm ] 



(49) 



where ( tr 7 2 ) = Yliji'Yij 7«j )• From this, it is easy to check that a cancellation leads to (m(k)) = 1 + C(7 4 
We can then expand eq. (|44p up to the first non-vanishing contributions. We obtain 



V c (k)= Vf\k) lm(k). 



( (m(k) 



1 



(trf 



1-(A + C) 
1)A) + (A) 

d Acey 



i 



Vf\k) 

d 



dvf\k) 
dh\k 



^ 1 d 2 vf\k) 

T v f) {k) d(\nky 



+ 



dink 

d 2 



+ 



(C 2 ) <? 



2 d{\nkf 



Vf\k) 



V 

2(T~~ ' 'din A; ' 2 N ' ' d{\nk) 2 I C 



(0) 



(50) 



This equation was also found in [10j 0- Neglecting tensor contributions, this reproduces the results of 
log-enhanced corrections to the power spectrum due to scalar fluctuations given in [TT|[T2]. 

Taking another point of view, we note that eq. (|50|) can be obtained by expanding £, the curvature 
perturbation in uniform-energy-density gauge, in terms of 6(f), the scalar field perturbation in flat gauge. Up 
to quadratic order in 5(j), the relevant gauge transformation is the one between £ and ( n = — (H/<fi)5<J) given 
in eq. (A. 8) of |61j . If we focus on terms that are leading order in slow-roll and neglect terms vanishing at 
superhorizon scales, all IR divergences arising from the expansion in (A. 8) of [61] are captured by our result. 
However, the complete IR correction requires the inclusion of term ~ <50 3 . This can be realized using the 
5N formalism, and it was shown in [11] that an appropriately modified version of this formalism correctly 
computes the scalar part of eq. ([50]) (see also Sect. [3] of the present paper). 

For a weakly scale-dependent power spectrum, the explicit values for ( Q 2 ) and ( tr j 2 ) were already 
given in eq. ([lip . Using the definitions of the spectral index of curvature perturbations rt£, its running 



and the tensor-to-scalar ratio r = N, 



2e, 



1 



dlnV 



(0) 
C 



d 2 lnV^ 0) 



dink " s d(lnk) 2 

the leading order correction to the power spectrum can be written as 



V c (k) 



Vf\k) 



1 + 



1 



(nc - l) 2 + a ? 



4r 



(n 



C 



V. 



(trf) 
8( C 2 ) 



^ 0) (A;) \n{kL)\ 



(51) 



(52) 



The agreement of eqs. ()50p and (|52[) with |10] is a non-trivial check for our approach. 



7 Note that in [TD] the result for corrections due to tensors is expressed in terms of the quantity (jq S ) = \ 53ij(7« Jij ) 
i{tr-y 2 } . Tensor contributions to inflationary observables, using different methods, have been also considered in [631174] . 



Note that the authors of [10] chose a different parameterization of the power spectrum, namely V 



(0) 

C 



fc"-( fc )- 



This leads 



to slightly different numerical factors. For instance, d 2 vi /d(\nk) 2 — [(n — f) 2 + 2a]Vi°^ in their parameterization. 
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We learn that long-wavelength modes provide log-enhanced contributions to the power-spectrum that 
are suppressed by second order slow-roll parameters, and by a factor of ■ The latter is determined by 

WMAP to be Vf ] ~ 2.3 x 10" 9 [75]. 

Let us stress the quite non-trivial fact that scalar and tensor long-wavelength modes contribute at the 
same (second) order in a slow-roll expansion. This is due to the cancellation leading to (m(k)) = 1 + 0(j 4 ). 
In any case, this property is specific of the power spectrum: as we will learn in the next section, corrections 
to non-Gaussianity parameters do not share this property. 

Having calculated the leading order correction to the power spectrum, we turn to evaluating scalar 

perturbations to all orders as described in section [2j We start with a inflationary potential for which the 

spectral index is constant. This is realized for the famous example of power law inflation [7B]. The potential 

/¥ 
q 

equations for scalar fluctuations can be solved exactly without having to rely on a slow-roll approximation. 
For this particular model, we assume that our scale of interest k is much smaller than /c max , reflecting the 
transition to scales that have never left the horizon during inflation. Hence, the integral in eq. (|18p is well 
approximated by setting £ m ; n to — oo. The power spectrum of curvature perturbations reads [771 175] 

/ , \ -2/(q-l) 



is V = Vq exp 



with constant q, and the scale factor evolves as a(t) = aot q . In this set-up, the 



pW(k) = K( k ^{k-J ■ ( 53 ) 
So the spectral index — 1 = —2/{q — 1) is constant as desired. We then obtain 

Vf\ke-Z) =vf\k) e - ( "c- 1)f " . (54) 
Plugging this expression into eq. (|18|) . one finds an integral that can be solved analytically. We get 

V ( (k) = Vf\k) exp ^ ( " C 2 ~ 1)2 ) • (55) 

This expression captures at all orders the contributions of long-wavelength modes. In a sense, we are 

providing the function whose series expansion has been found in |10j . Notice that the corrections are not 



independent of the scale k since is a function of k (e.g. = ln(fcL) for a weak scale-dependence of 
if). 

We expect a similar behavior including other contributions in more general models of inflation, for 
example associated with the running of the spectral index; however, solving the integral analytically might 
be more difficult in these models. For instance, the chaotic potential investigated in |10j . V(4>) = \<j) a 
(a > 0), leads to the following tree-level power spectrum: 

?(0)m _ ( N 4' H \ 2 _ 1 X A a+2, 



f'W = {-%-) {k) ■ (56) 

The scalar field value in dependence of the horizon-exit time of the mode k is given by the differential 
equation d(p/(dlnk) = —a/(j>. This can be integrated to yield 



^) = ^ 2 (A; max ) + 2aln^p . (57) 

Note that the condition C, > Cmin guarantees (j)(e"^k) > (f)(k mBX ) for all possible values (. Hence, the 
integral in eq. ([TH]) is well-defined and finite. As already described in sec. [21 the series expansion can be 
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recovered easily from the integral expression in eq. (|18p . For this purpose, one can expand vj- (e ^k), as 



given in eq. (j56[) . and make use of the moments for the Gaussian probability distribution (see eqs. (|16j) and 
(|17p ). Derivatives of the power spectrum (|56p w.r.t. ln/c can be expressed in terms of the spectral index 



n 



1 = —a(a + 2) /(j> and the model parameter a : 



d{\nk) 1 



+ 1 



(58) 



Hence, the series expansion is given by 



V c {k) 



Vf\k) 



(2n 



n=l 



(2ra)! 



^ 4 (J\ n 



(59) 



In the last part of this section, we discuss the question of convergence of the series expansion returning 
to the general case (see also [ID])- The series expansion in eq. (fTU|) was 



V c {k)=vf\k) 



(C 



2n 1 



d 2n Vf\k) 



n=l 



(2n)! d(lnfc) 2 " 



(60) 



and we have parametrically ( Q 2n ) ~ ( Q 2 ) n . Since In k! = In k — f, this is similar to a Taylor expansion of 
the power spectrum in In k around the scale k. At every order in the expansion, corrections consist of two 
counteracting contributions. On the one hand, there is the factor ( Q 2 ) n , with 



ce 



1/L 



(61) 



that involves a log-enhancement (even though it is suppressed by the smallness of the power spectrum). On 
the other hand, there are derivatives of the power spectrum that consist of slow-suppressed quantities and 
to which we will refer as late-time suppression. The word 'late-time' indicates that derivatives of the power 
spectrum in eq. ([60]) . i.e. the slow-suppressed quantities, are evaluated at the scale k. By contrast, the 
quantity ( C, 2 ) receives contributions from all modes in the range from 1 /L to the scale k. Convergence of 
the series expansion depends on the ability of the late-time suppression to compensate the log-enhancement 
due to (C 2 )- 

We will perform an order of magnitude analysis and, hence, we do not distinguish between quantities that 
are of the same order in slow-roll, like the slow-roll parameters e and rj. Instead, we generically characterize 
the slow-roll suppression by an appropriate power of e. A derivative d/d(lnk) acting on precisely 
corresponds to one such power of e. Hence, the late-time suppression is given by 



d 2n vf\k) 



v p {k) d{\nk) 



2n 



e 2n (k) 



(62) 



Therefore, in our order of magnitude analysis, eq. (|60p can be written as 



and the convergence of the series expansion requires e 2 (C 2 ) < 1 • 



(63) 
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Let us first consider a weakly scale-dependent power spectrum. Here, 'weakly scale-dependent' means 
that the power spectrum has only a very mild scale-dependence on the complete range from 1/L to k 
such that eq. f)61 1) essentially yields 

(C 2 }= Vf ln{kL) ~ — ln(fcL) . (64) 

The logarithm is given by the number of observed e-foldings N ~ Ht. Therefore, it remains to verify the 
relation 

e 2 (C 2 ) ~eHh < 1 . (65) 

As shown by [79 ,180j. the requirement of being in a non-eternally inflating phase constrains the time to obey 
t < RS ~ H~ 3 , where R and S are deSitter radius and entropy, respectively. Hence, the criterium for 
convergence reduces to e < 1, which is fulfilled by construction in slow-roll inflationary models. Therefore, 
under the assumption of a weakly scale-dependent power spectrum, the series is always converging. How- 
ever, this is not surprising. The assumption of a 'weakly scale-dependent' power spectrum can be made 
mathematically more precise by demanding that the scale-dependence of is negligible in the integral in 

eq. dSTj (such that ( ( 2 ) = 7^ 0) ln(fcL) ). This yields 

1 d n vf ] 

from which we could have concluded the convergence of the series in eq. (|6(jp directly. 

In spite of all that was said above, convergence breaks down for the model of chaotic inflation character- 
ized by the power spectrum (|56p (see also [10j ) - The integration in the expression of ( £ 2 ) in this model needs 
to be performed over several orders of magnitude in the scalar background field cp. Therefore the slow-roll 
parameter e ~ l/<^> 2 is changing over several orders of magnitude. This clearly violates the approximation 
of a weak scale-dependence. Consequently, convergence is not obvious in this model and an investigation 
of the convergence behavior requires a more precise evaluation of (C 2 )- Indeed, ( £ 2 ) ~ J(H 2 /e)dq/q is 
completely dominated by contributions at very early times ij. Hence, the expansion parameter ( Q 2 ) is much 
larger than vf* (k) In(kL). By contrast, the coefficients in eq. (|60p . i.e. the late-time suppression, consist of 
slow-roll parameters which are large compared to those at early times ti. Hence, the late-time suppression 
cannot compensate the enhancement at early times leading to a breakdown of convergence. 

In principle, this effect is also present for tensor corrections, though less severe since the power spectrum 
of tensor modes is not enhanced by 1/e. Therefore, in this model the breakdown of convergence due to 
scalar contributions occurs first. This observation is also in agreement with the findings in [ID], showing 
that the effect of scalars dominates. 

We note that a breakdown of convergence implies that one cannot trust conventional perturbation 
theory. However, this only applies to the conventionally defined power spectrum at sufficiently large L. In 
our philosophy, one should instead consider higher-order corrections to IR-safe quantities like the power 
spectrum V^\k) defined in eq. ([DJ. We know that the leading-order corrections to this object will not 
be log-enhanced. While we have not shown this in the present paper, we expect that also higher-order 
corrections will benefit from our IR-safe definition and hence that conventional QFT perturbation theory, 
based on the smallness of £ and of slow-roll parameters, will be reliable. 



[\n(kL)] n < 1 n >0 , (66) 

A' 



15 



5 Three-point function and the bispectrum 



The bispectrum accounts for the simplest contribution to non-Gaussianity. Starting from the three-point 
function in momentum space, one extracts the bispectrum from its connected part: 

(CfelCfe-Cfc-) = (27rf5^(k l + k 2 + k 3 ) B c (h, k 2 ) . (67) 

In this section, for definiteness we focus on non-Gaussianity of local form (see |81| for a recent review) Q 
The corresponding bispectrum is well-described by 

B c (h,k 2 ) = -fKL(h,k 2 ) [Pc(h)Pdk 2 )+ perms] , (68) 

where /nl is a slowly- varying function and we have introduced the uncurly power spectrum P^(k) = 
2tt 2 V^(k)/k 3 . Here and henceforth, we indicate with perms all non-trivial cyclic permutations of ki,k 2 
and k% = —(ki + k 2 ) . The dependence of P((k) on the long- wavelength background modes is characterized 

P c (k) = ^ e - 3 ?P c (0) (&')) . (69) 

With the formalism of sec. [2] for higher correlation functions, we may immediately write down /nl including 
long-wavelength corrections: 

, 5 B ( (k u k 2 ) 

JNL = ~ 



6 [P^(ki)P^(k 2 ) + perms] 
5 ( e~^ Bf\k[,k' 2 \ 



6 ( e- 3 C Pf ) {k' 1 )){ e- 3 C P^\k' 2 ) ) + perms 



(70) 



The remaining task is to evaluate (|70p . This requires knowledge on the tree- level bispectrum , which is 
model dependent. 

As an illustrative example, we consider the form 



B 



(0) 



6 



■ ^(k^P^ih) f ( (k 3 )+per7ns] . (71) 

This tree-level bispectrum is motivated by a curvature perturbation which is given by a Gaussian part, C G , 
plus fc_{k) times the Gaussian part squared, i.e. 

& = + fdk){C G H G )^ k ■ (72) 

Here, the operator * denotes a convolution. In concrete examples, depends on the scales k only by means 
of the dependence on times of horizon exit for each mode j82|184j. Note that the tree-level bispectrum 
(17111 has a slightly different scale-dependence than eq. (|68p . They only match for the popular assumption 
of being scale-invariant or in the squeezed limit, where one scale is much smaller than the others (say 
k\ <C k 2 ,k 3 ). Indeed, one has = /<; in these cases. 



We stress that the bispectrum ()7ip neglects the presence of intrinsic non-Gaussianity in the second 
order scalar field fluctuations. To include this contribution, one has to apply the bispectrum given by 
Maldacena [61]. In order to keep equations simple, we will neglect this presence of intrinsic non-Gaussianity 
and apply eq. ([7T]) . 



Other forms of non-Gaussianity can also be described with techniques similar to the ones we are are going to develop. 
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Note that the primary field of validity of eq. (|7ip is in multi-field, e.g. curvaton-type, models with 
observable non-Gaussianity. In addition, it arises in the squeezed limit of single field slow-roll inflation. 
In that case, our modified 5iV-formalism, presented in section [31 reproduces the correct result for fj^i i- e - 
/nl = 5/12 (1 — nf) (see appendix IB]) . This agreement shows that, contrary to the conventional 5N- 
formalism, our modified version of SN provides correct results also for the 3-point function in the squeezed 
limit. Consequently, the following calculation is correct in the squeezed limit, even though we made the 
simplifying assumption of negligible intrinsic non-Gaussianity. 

Proceeding as we did for the power spectrum, we perform a slow-roll expansion for the quantities inside 
the averages in eq. (|70p . focussing on the non- vanishing contributions at leading order in slow-roll. After 
some calculation, this yields 



f _ AO) 

JNL - J NL 



AO) 

JNL 



AO) 

JNL 



1 + 



Qjh, fc 2 , fc 3 ) PP(h) PP(k 2 ) f c (k 3 ) + perms 
P f (0) (h) P c (0) (k 2 )f c (k 3 ) + perms 



k 2 , h) A:f + 0(fc 3 , fci, k 2 ) k'l + fi(fc 2 , fc 3 , fcx) kf 



j_ h.3 _i_ L.3 



Bf\k 1: k 2 ) 



6 P f (0) (fci)P c (0) (k 2 ) + perms 



20 



(trf ) 3(A;i • k 2 ) 2 - 1 



-(trf)<(2 3(fe-AaV 



1 dP, 



(0) 



(73) 
(74) 
(75) 



-p(o) din A; 
1 C 



+ 3 



^ /cd(lnA;) 2 



(o)' 



din A 



(h ■ k 3 ) 2 + (fe • k 3 ) 2 - l 



(o) 



f^d\ak 



+ 4 



1 d/ c g 
din A; din A; 



(76) 



Here fej • A;.,- corresponds to the cosine of the angle between the vectors k{ and kj . ( tr 7 2 ) and ( £ 2 ) are defined 
as before. In the previous expression for Q, the scale at which we evaluate fc_ and their derivatives is 
any one of the k\ : the difference among these quantities evaluated at different scale is slow-roll suppressed 
with respect to the contributions we are examining. As a result, a cancellation of corrections originating from 
the numerator and the denominator in eq. (|73p occurs. This removes terms containing second derivatives 
of the power spectrum. Moreover, it is sufficient to take into account the naive scaling P^°\k) ~ A;" 3 in 
eq. (|73p . leading to the simpler form in eq. (|74p . In eq. (|75p . we defined the leading order non-Gaussianity 
parameter 

In some cases, it may be useful to perform an average over directions of the vectors. However, we note 
that the <5-function sets constraints on this averaging procedure. As an example, we focus on the particular 
case of squeezed configurations, i.e. k\ <C k 2 ,k 3 . For these configurations, one of the permutation terms can 
be dropped and the 5-function requires {k 2 ■ k 3 ) 2 = 1. The pair of unit- vectors k\, k 2 or k\, k 3 is statistically 
independent. Hence, the directional averaging gives (ki-k 2 ) 2 = (k\-k 3 ) 2 = 1/3. Therefore, having performed 
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the directional averaging, the expression for squeezed configurations reads 



f _ AO) 

JNL - / NL 



1 + Q(k\, k 2 , k 3 



+ 



20 V ' f^dlnk 

(C 2 ) / i ^ 2 /c 



f<d(lnky 



+ 2 



<p(o) dink 
1 C 



+ 4 



dV, 



(0) 



fc d In k dink 



(77) 



(78) 



We stress that under the directional averaging the first term on the right-hand side in eq. (|76p vanishes. 
Indeed, we will see below that, keeping the directional information, precisely this term turns out to be the 
leading order correction. Therefore, this example illustrates that such procedures have to be handled with 
care. 

Neglecting tensor fluctuations, the special case of corrections to /nl hi squeezed configurations was also 
discussed in [10] , In this case, corrections are solely given by the last line in eq. (|76p . which reads 



n(ki,k 2 ,k 3 



(C 2 ) 1 dPf ( 



/Cd(lnfc) 2 



+ 2 



1 dV 



(o) 
C 



1 C 



+ 4 



1 



df C 



dV, 



(o) 
C 



j?^ <p(p) dink dink 



(79) 



Our result basically agrees with the findings of [10]. In order to have complete agreement, one needs to take 
into account the runnings of and in the calculation of the bispectrum in [10] ( their eqs. (5.9)-(5.11) ). 
Including this effect, their final formula for /nLj eq. (5.14), slightly changes. Like in our findings, the running 
of the power spectra from the denominator disappears since it cancels against corresponding terms from 
the numerator. The effect of the running of fc_ appears precisely as the first term on the right-hand side 
of eq. (|79p . The second and third term on the right-hand side of eq. (|79|) are already present in eq. (5.14) 

in nm. 



We now return to the general form of corrections to /nl, i-e. eqs. (|73p - (|76p . Remarkably, we find that in 
single-field, slow-roll inflation tensors provide the dominant contribution in slow-roll, i.e. the first term on 
the RHS in eq. (|76p . This contribution results in a correction proportional to first order slow-roll parameters, 
while the others are of second order. Indeed, we observe that, at leading order in slow-roll, the dominant 
contribution is originating from the prefactor 



{m(k 1 )m(k 2 )) 



1 



(80) 



which multiplies the tree-level bispectrum. Prom this we find that the dominant log-enhanced contribution 
to /nl, in a slow-roll expansion, reads 



f _ AO) 
/nl - / NL 



1 + £ 7f } In(fcL) (3(fcl • h) 
5 ^ 



1) fcf + (3(fc 3 • ^i) 2 - 1) kl + (3(fc 2 • k? - 1) k 



fe3 i US 
1 2 



1.3 
ft 3 



(81) 



Interestingly, these log-enhanced contributions to /nl do not depend on the tilt of the power spectrum, and 
so are also present for spectral index equal to one. 

In conclusion, log-enhanced contributions to /nl can be expressed in terms of observable quantities. 
Tensor contributions are proportional to first order slow-roll parameters, and are suppressed by the tree- 
level power spectrum. Very similar results hold for parameters associated to the trispectrum, $nl an d tnl- 
It is straightforward to obtain them proceeding exactly as done in this section. 

Local non-Gaussianity in single field, slow-roll inflation turns out to be small. On the other hand, models, 
in which a second field takes part in the generation of curvature perturbations as in the curvaton scenario, 
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can lead to large values of /nl (see e.g. [85] )• In the approximation in which only the curvaton field is 
responsible for curvature perturbations, the tree level bispectrum reads 



fa) = fa(h,fa) 



~~kj 



2tt 



V a {ki) -pr-Vaifa) + perms 



(82) 



where a indicates the curvaton field. Our approach can be applied also to this case, although it requires 
additional work to calculate the contributions of long-wavelength scalar modes to inflationary observables, 
since more than one scalar field is present. We outline a method to do this in Appendix [Bj but a more 
complete discussion of this issue is left for future work. In this case, enhancement effects associated with 
long-wavelength modes could turn out to be more important than the ones discussed so far. 



6 Conclusions 

We have considered IR effects associated with backreaction of long-wavelength scalar and tensor modes 
in inflationary backgrounds. We proposed an infrared-safe definition of correlation functions involving 
curvature fluctuations, with no sensitivity on long- wavelength contributions. The essential idea was to make 
use of the proper invariant distance on the reheating surface where the curvature perturbation is evaluated. 
By using the invariant distance, one automatically absorbs longer wavelength modes in the background and 
obtains ra-point functions for the curvature perturbation that are free from IR contributions associated with 
long-wavelength modes. We showed how to re-interpret our results in terms of conventionally defined n-point 
functions. This allowed us to provide closed expressions for the latter that manifestly exhibit the dependence 
on long-wavelength modes. In our approach, IR corrections automatically emerge in a resummed, all-orders 
form. We then applied our approach to the analysis of inflationary observables built from (conventionally 
defined) two- and three-point functions of the curvature perturbation. We showed how to compute the 
leading scalar and tensor IR effects on the power spectrum and on the bispectrum, in single field, slow-roll 
inflation. Our corrections to the power spectrum (both from long-wavelength scalar and tensor modes) 
and to /nl (from long-wavelength scalar modes) agree (essentially) with Giddings and Sloth [10] (obtained 
by somewhat different methods). The advantage of our approach is that it directly provides resummed, 
all-orders expressions. We extend [TU] by tensor corrections to /nl- This is, in fact, the dominant piece! We 
also explicitly computed, in a specific inflationary model, the complete, all-orders expression for scalar long- 
wavelength contributions to inflationary observables. Furthermore, we analyzed the question of convergence 
of IR corrections. Using entropy bounds given in [79l[80], we found that for a weak scale-dependence the 
convergence of the series of IR corrections is guaranteed. However, despite the existence of these entropy 
bounds and the fulfillment of slow-roll conditions, the convergence of the IR-correction series may break 
down if the scale-dependence is not sufficiently weak. 

Summarizing, we have provided a simple formalism to calculate and investigate inflationary IR correc- 
tions. Maybe more importantly, we have provided simple definitions of IR-safe correlation functions which 
make it possible to avoid IR enhancement altogether. 

We have also shown that in all cases, where the <5iV-formalism is applicable, our results can be equivalently 
obtained in terms of a suitable generalization of the 5N- formalism, extending the discussion of [TTJ. In the 
present work, we included the effects of graviton long-wavelength modes, and we explained how to calculate 
IR contributions to arbitrary n-point functions involving curvature perturbations. 

A natural question is how to extend our results to the case in which more than one field plays an active 
role in generating the curvature perturbations. In this case, IR effects might play a role more important 
than the one for single field inflation. We outlined in an Appendix a method to treat this problem, but we 
leave a more complete discussion for future work. 
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A Comparison of and V ( 



(0) 

C 



The definitions of the IR-sensitive power spectrum and the IR-safe power spectrum yf are 
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In the first line, we included the integral and prefactors in the average. Note that this does not affect the 
averaging process over pairs of points separated by the coordinate vector y. From the second to the third 
line, we performed a coordinate transformation of the integration variable from y to z. Since the determinant 
of e 1 is one, tensor fluctuations do not effect this transformation. Therefore, only scalar fluctuations appear 
as a prefactor in the third line. Consequently, we need to add tensor fluctuations by hand in this prefactor, 
in order to express the third line in terms of the IR-safe power spectrum . This results in the prefactor 

in the last line, which only consists of tensor fluctuations. In this last line, the vector A; is a unit vector in 
/c-direction and the average is performed over the background quantities £(x) and jij(x). 



B Extension of £7V-formalism 

In this appendix, we discuss in more detail how our results can be understood in terms of a 5N approach. In 
a previous paper written in collaboration with Byrnes and Nurmi, we showed how a suitable extension 
of the (WV-formalism allows for the computation of leadingdog contributions to the power spectrum, due 
to scalar long-wavelength fluctuations. Here, we extend our work to include tensor modes and to compute 
log-enhanced corrections to non-Gaussianity parameters. 

Let us start with single field inflation. For this purpose, we will adopt the same gauge as in section [3) 
By means of the 5iV-formalism, the curvature perturbation £ can then be expressed in terms of the number 
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of e-foldings evaluated on a background given by the scalar field <j> and its perturbation 



Q = N [cf> + 5cf>] - (N) . 
The previous expression admits an expansion in terms of scalar fluctuations 

= N4k')5^(k',g ab ) + \N H {k') [{50*50) % {k',g ab ) - (50*50)] + 



(89) 



(90) 



where we use the notation of the main text, and the Gaussian, first-order scalar fluctuation 5^0 is given in 
eq. (|34p . The effect of long- wavelength modes is encoded in the shift of the time of horizon exit, and in the 
function m, contained in the expression of 5r(f). Long-wavelength mode contributions are controlled by the 
averaged quantities £ and 7^-. 

With the previous expression, we neglect intrinsic non-Gaussianity of 5^0. Hence, this formalism is only 
applicable in situations where this is negligible. However, this condition is fulfilled in several models, e.g. 
models in which non-Gaussianity of the local form can acquire sizeable values as in multiple field inflation or 
curvaton-like mechanisms [S3]. In light of these models, it is worthwhile to develop formalisms that neglect 
the presence of second order fluctuations. 

Using the previous formula, it is straightforward to compute n-point functions of curvature perturbations, 
and compute leading log-enhanced corrections to inflationary observables. It is important to stress that, due 
to the choice of a gauge with 5^j) = for q <C k, convolutions appearing in the second term of eq. ([90]) do 
not involve integration over all the modes, but have a lower cut-off slightly below the scale k. This implies 
that convolutions, when appearing in n-point functions, do not provide further log-enhanced contributions 
with respect to the ones associated with long-wavelength background modes. All the IR dependence is then 
contained in the quantities £ and 7. 



As an example, let us work out explicitly the expression for the three-point function in the squeezed 
limit, including the effects of long- wavelength modes, using eq. (|90p . Our method is similar to [SS]. The 
second term on the right-hand side in eq. ()90p is irrelevant for squeezed configurations and will be neglected. 
The contribution to the bispectrum is 

<C^ 2 Cfc 3 > = ^=p ([m3(fc 1 )(iVr)(fc , 1 )aJ [m3(fe2)(iV^)(^)aJ [m^k) (N^H)(k^) a J ) . 

(91) 

In the limit in which k\ <C k 2 ~ k^, the size of the vector k\ is comparable to the size of the long- wavelength 
modes relative to the vectors k 2 and k 3 . The latter are included in the shift of the time of horizon exit ty 

and on m^kz), respectively, t^/ and m(k^). Taking into account this fact, and using Wick's theorem, we can 
write in this limit the following non-vanishing contribution 



(C^A i ) = (^ {3 Hh + k 3 ) 
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(92) 



Since k\ <C k 2 , the only possibility to contract is the background contribution originating from (N < f > H) 2 (k' 2 ) 
By expanding the latter and by means of the definition of this contraction yields 
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Therefore, we find for the bispectrum: 
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At leading order, we neglect the contribution from the background (k^ — > ki and m(ki) = 1). This provides 
the following tree-level result for the non-Gaussianity parameter 



L ,-> (0) 



5 dlnV K ( _ 

k 



12 dink 



^(1-«C(*)) > (95) 

that is Maldacena's consistency relation |61|I62|. Hence, the complete form for /nl> obtained from formula 
(fM|) . can be expressed through f^, giving 

(mfavjPW) m{h)vf\k> 2 ) 

JNL = ; 77Ts ; 77a , (9b) 

in agreement with formula (|70p in the squeezed limit. We can then proceed as done in the main text to 
extract leading log-enhanced contributions. 

Let us briefly discuss the case in which multiple scalar fields affect the curvature perturbation. The 5N- 
formalism is very well suited to study this case, as discussed in the original paper by Sasaki and Stewart [72] . 
We adopt a gauge with vanishing scalar metric fluctuations, i.e. the long-wavelength modes of the scalar field 
are not vanishing. The curvature perturbation (that in the case of multiple fields is not generally conserved) 
can be expressed as an expansion in terms of all the scalar fields involved 

C n (t f ) = N![t f ,{M} + {(50 7 *50 J ),:[{0o}]- (<50 J *<50 J )} + --- , (97) 

where the capital latin indices of N denote derivatives w.r.t. the scalar fields and summation over repeated 
indices is understood. Here, {0o} denotes the dependence on the homogeneous values 0q of the scalar fields. 
As in Sasaki and Stewart, we have replaced the dependence on the time of horizon exit tf- of the various 
functions, with the value of homogeneous solutions of the scalar equations at ty. 0n = 0o(^fc) • 

Then, the inclusion of the effects of scalar and tensor long-wavelength modes can be done as in the 
previous sections, although the procedure is a bit more laborious. We express a given function 7 (£, x), the 
solution of the field equations, as 

7 (i, x) = 4(t) + <50 J (i) + 5(f) 1 (t, x) , (98) 

where 5(f) 1 (t) is an average over long- wavelength modes, similar to the ones we performed in the main text. 
The effect of long- wavelength scalar fluctuations is to shift the values of 0q, that appear in eq. ([97]) . to 
4>q + 8<j> . In a sense, they play the same role of shifting the time of horizon exit, although with multiple 
fields there is not a one to one correspondence between time and values of the scalar solution. After passing 
to momentum space, the inclusion of long-wavelength scalar perturbations implies that the expansion in 
eq. (f§T|) becomes 

c n (t f ) = Nj [t f , {4> + m] b^ 1 [{&, + h}} 

+ i N u [t f , {0o + 50}] { {H 1 * H J ) £ [{0o + 50}] - (50 J * <50 J ) } + • • • (99) 

to take into account the shifts of the homogeneous solution of the scalar fields. The inclusion of tensor long- 
wavelength contribution, on the other hand, is very simple: since correlations between tensor and scalar 
modes vanish, the effect of tensors is precisely identical to that discussed in the previous sections. It can be 
taken into account with a proper redefinition of the time of horizon exit of a given mode tf. —> ty. One can 
then repeat in this case the very same steps that we took in the previous sections, generalizing our results 
to multiple fields. This will be done in future work, where we will also discuss in this context the possibility 
of having large non-Gaussianity from loop effects [33], with sizeable scale-dependence 
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